function pi = pi_f(pm,x,e,p,w,t,r)

pi = (1-pm.rho).*(1-t).*p.*(x.*e).*l_f(pm,x,e,p,w,t,r).^pm.rho;